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Recent Advances in Reaction-Diffusion Equations 
with Non-Ideal Relays 

Mark Curran 1 2 3 4 5 Pavel Gurevich 2,3 Sergey Tikhomirov 4,5 


Abstract 

We survey recent results on reaction-diffusion equations with discon¬ 
tinuous hysteretic nonlinearities. We connect these equations with free 
boundary problems and introduce a related notion of spatial transver- 
sality for initial data and solutions. We assert that the equation with 
transverse initial data possesses a unique solution, which remains trans¬ 
verse for some time, and also describe its regularity. At a moment when 
the solution becomes nontransverse, we discretize the spatial variable and 
analyze the resulting lattice dynamical system with hysteresis. In par¬ 
ticular, we discuss a new pattern formation mechanism — rattling , which 
indicates how one should reset the continuous model to make it well posed. 


1 Introduction 

1.1 Motivation 

In this chapter we will survey recent results on reaction-diffusion equations with 
a hysteretic discontinuity defined at every spatial point. We also refer to m 
and the more recent surveys by Visintin [4ll5| for other types of partial differential 
equations with hysteresis. 

The equations we are dealing with in the present chapter were introduced 
in HE! to describe growth patterns in colonies of bacteria (Salmonalla typhir- 
murium). In these experiments, bacteria (non-difussing) are fixed to the surface 
of a petri dish, and their growth rate responds to changes in the relative concen¬ 
trations of available nutrient and a growth-inhibiting by-product. The model 
asserts that at a location where there is a sufficiently high amount of nutrient 
relative to by-product, the bacteria will grow. This growth will continue until 
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the production of by-product and diffusion of the nutrient lowers this ratio be¬ 
low a lower threshold, causing growth to stop. Growth will not resume until the 
diffusion of by-product raises the relative concentrations above an upper thresh¬ 
old that is distinct from the lower. Numerics in IB] reproduced the formation 
of distinctive concentric rings observed in experiments, however the question of 
the existence and uniqueness of solutions, as well as a thorough explanation of 
the mechanism of pattern formation, remained open. 

Another application in developmental biology can be found, e.g., in [8], and 
an analysis of the corresponding stationary solutions in [5]. 


1.2 Setting of the Problem 

In this chapter we will treat the following prototype problem: 


ut = Au +f(u,v), v = n{^ 0 ,u), 
u|t=o = ¥>, x G Q, 


du 

dv 


= 0 . 

O'Qt 


{x,t) G Qt, 


(1) 

( 2 ) 

( 3 ) 


Here Q C R" is a domain with smooth boundary, Qt := Q x (0, T), where 
T > 0, O'Qt '■= dQ x (0 ,T), u is a real-valued function on Qt , and %(£o>u) is 
a hysteresis operator defined as follows (see Figure [Ur). Fix two real numbers 
a < /3, an integer £ 0 G {—1,1}, and two continuous functions H\ : (—oo, j3\ —> R 
and H_i : [a, oo) —> R such that Hi (it) ^ H_i(it) for u G [a, j3\. Define the sets 

Ei := {{u,v) G R 2 | u G (— oo,/3),v = 

E_i := {(it, it) Gl 2 | it G (a,oo),it = H_i(it)}. 


Definition 1.1 Letu,v : [0,T] —> R, where u is a continuous function. We say 
that v = H(<!;o ; u ) tf the following hold: 

(1) (u(t),v(t)) G Ei U E_i for every t G [0,T]. 

(2) //it(0) G (a,/3), then v(0) = H^ o (u(0)). 

(3) If it(fo) G (a, /3), then v(t) is continuous in a neighorhood ofto- 

The operator 0 ;^) is called the non-ideal relay and item [3] means that the 
non-ideal relay jumps up (or down) when it = a (or u = /3). This defini¬ 
tion is equivalent to the definitions of non-ideal relay found in mmm- k 
"H(£ o, u)(t) = Hj(u(t)), then we call £(£) := j the configuration of H at the 
moment t , and we call £o the initial configuration. Now let it : Qt —> R be 
a function of [x,t) and '■ Q {—1,1} a function of x, then T-L{f,o,u)(x,t) 
is defined in the same way by treating £ as a parameter, i.e., there is a non¬ 
ideal relay at every x G Q with input u(x,t), configuration £(x,t), and initial 
configuration fo{x). 
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Figure 1: The hysteresis operator with H\{u) = 1 and H-\(u) = — 1. 


1.3 Set-Valued Hysteresis 

First results on the well-posedness of © © were obtained in mini for set¬ 
valued hysteresis, and their model problems are worth explaining in more detail. 
In both papers, the uniqueness of solutions as well as their continuous depen¬ 
dence on initial data remained open. 

First we discuss the work of Visintin |j]T, which treats © © for arbitrary 
n > 1 with replaced by a set-valued operator called a completed relay 

(see Figure©}). We still use the thresholds a < /3, and will consider constant 
hysteresis branches H\{u) = 1, and H_i(it) = —1. We also define the set 
E 0 := {( u,v ) £ R 2 | u £ [a,ff\,v £ (—1,1)}. 

Definition 1.2 Let u,v : [0,T] —> R, where u is a continuous function, and let 
£o £ [—1,1]. We say v £ 'Hvis(£o,w) if the following hold: 

(1) ( u(t ), v(t)) £ Ei U £_i U So for every t £ [0, T]. 

(2) If u( 0) £ (a,/3), then u(0) = £o! if u(0) = a (or ft), then r>(0) £ [£o, 1] (or 
r(0)£[-Uo]). 

(3) If u(to) £ (a,/3), then v(t) is constant in a neighborhood of to. 

(4) If u(to) = a (or /3), then v{t ) is non-decreasing (or non-increasing) in a 
neighborhood of to. 

By treating i as a parameter, 'Hvis(£o,u) is defined for u : Qt —> R as we 
have done previously for 7-L(fo,u). Visintin [T3] proved the existence of u and v 
such that the equation 


u t = Au + v, v £ 'Hvis(£o,M), 

with n > 1, Dirichlet boundary conditions, and initial data tp is satisfied in a 
weak sense in Qt- Visintin [T3] and more recently Aiki and Kopfova [141 proved 
the existence of solutions to modified versions of ia ei, where the hysteretic 
discontinuity was a completed relay responding to a scalar input. A non-ideal 
relay with vector input, as in El El, behaves almost identically to a non-ideal 
relay with scalar input, but for clarity of exposition we only consider scalar 
inputs in this chapter. 
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Let us now turn to the model hysteresis operator %Ait(Co,u) proposed by 
Alt in [12] (see FigureflJ). We still consider Hi{u) = 1 and H_\(u) = —1, and 
introduce the set 

Eq : = {(«, v) e R 2 | u = a, v G [—1,1)} U {(«, v) £ ! 2 j u = ft,v € (—1,1]}. 


Definition 1.3 Let u,v : [0,T] — > R, where u is a continuous function, and let 
Co £ {—1,1}. We say that v £ HAit(fo,u) if the following hold: 

(1) ( u(t), v(t)) £ Ei U E_i U Eo for every t £ [0, T]. 

(2) Ifu{ 0) £ [a, ft], then v(0) = Co- 

(3) If u(to) £ (a,/3), then v(t) is constant in a neighborhood of to. 

(4) If u{to) = a (or ft), then ift) is non-decreasing (or non-increasing) in a 
neighborhood ofto- 

One can define %Ait(£o,w) for u : Qt —> M by treating x as a parameter as we 
did when defining o,u) and Ffvis(Co,w). 

To highlight the main difference between the completed relay TIvis (Co , u) and 
Alt’s relay ^Ait(Co,u), suppose that ^ V is(Co,w)(fo),^Ait(Co,u)(fo) £ (-1,1) 
and u(to) = /3 has a local maximum at time to- Then, as soon as u decreases, 
Hah jumps to —1, however Hy is remains constant. 

Let us introduce the notation {u = a} := {(x,t) £ Qt | u(x,t) = a}, 
with {u = ft} defined analogously. Alt’s existence theorem can, omitting the 
technical assumptions, be stated in the following way. Let n = 1 and suppose 
(ip, Co) £ Si U E_i. Then the following holds: 

(1) There exists u and v such that v £ 'HAit(Co, u ) a - e - dr Qt and 
u t = u xx + v a.e. on {(x, t) £ Q T \ u{x, t) £ {a, /3}}. 


(2) We have 


Ut — 'U'xx on {(x,t) £ Q t \ u(x,t) £ {a,/3}}, 
v £ [—1,0] on {u = ft}, and v £ [0,1] on {u = a}. 

(3) Items [5] 0 of Definition 11.31 hold in the following weak sense: 

For every ip £ C^°(Q x [0, T)) with ip > 0 on {Q x [0, T)} D {u = a} and 
ip < 0 on {Q x [0, T)} (~l {u = ft}, 


’ Qt 


(■ r) — iio)ipt dxdt < 0. 
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1.4 Slow-Fast Approximation 

Equations of the type m © are deeply connected with slow-fast systems where 
the variable v is replaced by a fast bistable ordinary differential equation with 
a small parameter S > 0 

5v t = g(u,v). (4) 

A typical example are the FitzHugh-Nagumo equations, where g(u , v ) = u + v — 
i> 3 /3 and the hysteresis branches H\(u) and 7T_i(it) are the stable parts of the 
nullcline of g (see Figure [2]). The question of whether the hysteresis operator 



Figure 2: a) The nullcline of the S'-shaped nonlinearity g(u,v). b) Hysteresis 
with nonconstant branches H\{u ) and H-\(u). 

approximates the fast variable v as S —> 0 has been addressed for systems of 
ordinary differential equations (see, e.g., [T51ITC] and further references in [IT]), 
however the corresponding question for partial differential equations is still open. 

1.5 Free Boundary Approach 

Problem (JTJ) |3]) with hysteresis has two distinct phases and a switching mech¬ 
anism, hence it can be considered as a free boundary problem. First observe 
that the hysteresis H naturally segregates the domain into two subdomains 
depending on the value of £(x,t). Denote 

Qj :={x€Q | £ 0 (x) = j}, j = ±1. (5) 

Let us look at how the free boundary Q\ fl Q-\ can evolve for a simple example 
on the interval Q = (0,1). Consider a neighborhood U of x £ Q, and suppose 
at time t = 0, Q\ fl U and Q -1 fl U are subintervals separated by a point b £ U 
(Figure [3ja) - Let u(x,to) > /3 for x < b, u(x,to ) < /3 for x > 6, and let x = a{t) 
be the unique solution of u(x , t) = /? in U. If at time t\ > 0 the value of u 
at points x > b have already risen above /3, then £(x,t) has switched from 1 
to —1. These are the points x such that b < x < a(fi) (Figure [3b). Now if 
at time ^2 > ti the value of u at the switched points has fallen below j3 again, 
i[x : t) remains switched. These are the points x such that a{t 2 ) < x < a(t \) 
(Figure (3t). More succinctly, (;(x,t ) = — 1 if x < b(t) and £(x,t) = 1 if x > b(t), 
where b(t) = maxo< s <t a(s). 

The point of this example is to illustrate that the free boundary does not 
in general coincide with the points where u is equal to one of the threshold 
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Figure 3: An example of the hysteresis configuration £ responding to an input u. 


values. This is different from the two-phase parabolic obstacle problem (see, 
e.g., mm), which © © reduces to if a = /?. 

Assume the derivative ifi'(b) in the above example was non-vanishing on the 
boundary {6} = Q\ fl Q~\. This is an example of transverse initial data, and 
whether the initial data is transverse or not will play an important role in the 
analysis of problem © ©■ 

1.6 Overview 

This chapter is organized in the following way. 

In Section [2] we will investigate the well-posedness of © © for transverse 
initial data. For n = 1 the existence of solutions and their continuous depen¬ 
dence on initial data was established in mi, uniqueness of the solution in [20] 
and the analogous results for systems of equations in [21]. Preliminary results 
for n > 2 were obtained in [22] . 

In Section[3]we consider the regularity of solutions u, in particular, whether 
the generalized derivatives u XiXj and Ut are uniformly bounded. We will sum¬ 
marize the results of [23], where the authors proved that these derivatives are 
locally bounded in a neighborhood of a point not on the free boundary. They 
also showed that this bound depends on the parabolic distance to the parts of 
the free boundary that do not contain the sets (w = a} or (u = /3}. 

In Section [4] we consider non-transverse data and the results of [24]. We 
will analyze a spatio-temporal pattern (called rattling) arising after spatial dis¬ 
cretization of the reaction-diffusion equation and discuss its connection with the 
continuous model with hysteresis operators "H, TLyis, and TLah- 
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2 Transverse Initial Data 


2.1 Setting of a Model Problem 

In this section we will discuss the well-posedness of problem © © under the 
assumption that <p is transverse with respect to £o, a notion which we will make 
precise shortly. In order to illustrate the main ideas, we will treat the following 
model problem in detail and then discuss generalizations at the end of this 
section (see Subsection m - Let /i_i < 0 < hi be two constants, and let the 
hysteresis branches be given by H\{u) = h± and H_i(u) = h- 1 . Consider the 
prototype problem 


u t = Au + H{f 0 ,u), 
w|t=o = <p, x £ Q, 


du 

dv 


= 0 . 

O'Qt 


(■ x,t ) e Q t , 


( 6 ) 

(7) 

( 8 ) 


We will treat n = 1 in Subsection 12.21 (see PU12Q3) and n > 2 (see [22] ) in 
Subsection 12.31 Throughout this subsection we will always assume that <p and 
£o are consistent with each other, i.e., if <p{x) < a (or ip(x) > (3), then £o(20 = 1 
(or £ 0 (a;) = —1)- In particular, this means that for every x £ Q, £,(x,t) is 
continuous from the right as a function of t £ [0, T). 

Since in general u) £ Lq^Qr), we will look for solutions in the Sobolev 
space Wq^^Qr) with q > n + 2. This is the space consisting of functions 
with two weak spatial derivatives and one weak time derivative from Lg(Qr) 
(see [23 Chapter 1 ]). If u £ Wq' 1 (Q T ), then for every t £ [0,T] the trace is 

well defined and u(-,t) £ Wq~ 2 ^ q (Q) (see, e.g., [23 P- 70]). To ensure that (p 
is regular enough to define the spatial transversality property, we henceforth 
fix a 7 such that 0 < 7 < 1 — (n + 2)/q. It follows that if ip £ Wq 2 ^ q (Q), 
then ip £ CF(Q) and Vip £ (C 7 (Q)) n , where C 7 is the standard Holder space 
(see [56] Section 4.6.1]). 

The subspace W 2 ^^ q (Q) C Wq 2 ^ q (Q) of functions with homogeneous Neu¬ 
mann boundary conditions is a well-defined subspace, and in this section we 
always assume that (p £ W 2 '~N /q (Q)- 

Definition 2.1 A solution to problem © © on the time interval [0,T) is a 
function u £ Wq :1 (Qr) such that ([S]) is satisfied in LqlQx) and u satisfies (0 
and m in terms of traces. A solution on [0, 00 ) is a function « : Qx [0, 00 ) —> R 
such that for any T > 0, u\q t is a solution in the sense just described. 


We note that if u £ Wq^^Qr), then r H((; 0 ,u) is a measurable function on Qt 
(see [T] Section VI. 1]). 


2.2 Case n = 1 

Let Q = (0,1) and Qj be given by ©. 
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Definition 2.2 Let ip £ C' 1 (Q). We say p is transverse with respect to £o if 
the following hold: 

(1) There is a b £ (0,1) such that Q—i = {x | 0 < x < 6} and Q± = {x \ b < 
x < 1}. 

(2) If p(b) = /?, then < p'(b) < 0. 

An example of <p and £o satisfying Definition 12.21 is given in Figure©,. 

Definition 2.3 A solution u is called transverse if for all t £ [0,T], u(-,t) is 
transverse with respect to £(-,i). 

Theorem 2.4 (See [11, Theorems 2.16 and 2.17]) Let p £ be 

transverse with respect to £o. Then there is a T > 0 such that the following 
hold: 

(1) Any solution u £ of problem © -© is transverse. 

(2) There is at least one transverse solution u £ W^ ,1 (Qt) of problem (fU jl ~(f% |l . 

(3) If u £ W^ ,1 (Q t ) is a transverse solution of problem (O ©. then it can 
be continued to a maximal interval of transverse existence [0,T max ), i.e., 
u(x,T max ) is not transverse or T max = oo. 

We will sketch the proof of Theorem 12.41 part [2] assuming that tp(b) = /3 and 
p'(b) < 0. 

Let us define the closed, convex, bounded subset of C[0,T] 

B := {b £ C[0, T) | b(t) £ [0,1], 6(0) = 6}. 


For any € B , define the function 


F(x,t) 


h- 1 if 0 < x < 6o(t), 
h\ if 6o(t) < x < 1. 


(9) 


Let u £ Wg' 1 (Qx) be the solution to problem © © with nonlinearity F in 
place of "H(£ o,u). We claim that T can be chosen small enough such that the 
configuration £,{x,t) of T-L(f;o,u) is defined by a unique discontinuity point b(t). 
Note that we do not yet claim that F = H(£q,u). 

To prove the claim, first fix To > 0. ft is a result of classical parabolic 
theory [25J Chapter 4] that for all T £ [0,T o ] 

IMIct(qT) 11 Ux 11 C~i(Qt) — (||^IU,(q t ) + — C 2 , ( 10 ) 

where Ci, C2,... > 0 depend only on Tq and q. The claim now follows from m 
with the help of the implicit function theorem. 

Observe that u is a solution of problem © © if r H(^o,u) = F, i.e., b 0 = b. 
We therefore look for a fixed point of the map 1Z : B —> B, F(b 0 ) := b. 


Consider b 0 i,b 0 2 £ B and define F\,F 2 via boi,bo 2 similarly to ©, and 
let u±,U 2 be the corresponding solutions. Observe that F\ ^ F 2 only if 
mm(boi(t),bo 2 (t)) <x < max(boi{t),bo 2 (t)), in particular, 

\\ui- u 2\\c-r(Q^) + \\uix-U2x\\ci(Q^) < Cl ||-Fl - F 2 ||i,(Q T ), 

< C3II60I - ^02||c[0,T]- 

Applying m again, and using <p'(b) > 0 and the implicit function theorem, we 
see that the left hand side of CD bounds ||ai — cl 2 \\c[o,t}- O ne can additionally 
show that ||ai — ^ 21 |c[o,t] bounds ||&i — 62 llofo.T ’]5 hence 

ll&l “ Mlc[ 0 ,T] < ll a l _ a 2 ||c[ 0 ,T] < C 4 II 601 - &02||,T]‘ ( 12 ) 

In particular (1121) shows that TZ is a continuous map on B. Moreover, one can use 
m to show that 7Z(B) is bounded in C 7 [0,T], and since C 7 [0,T] is compactly 
embedded into C[0,T], the Scliauder fixed point theorem implies that 1Z has a 
fixed point. 

Theorem 2.5 (see (2131 . Theorem 2.2]) If u 1 and, U2 are transverse solutions 
of problem ©-© with the same ip, then u± = U 2 - 

We prove the theorem by expressing solutions as a convolution with the Green 
function G(x, y , t, s) for the heat equation with Neumann boundary conditions. 
Let us use this function to estimate the solution w = u± — U 2 of the heat equation 
with zero initial data, Neumann boundary conditions, and the right hand side 
h = — H(fo,U2)'- 


\w{x,t)\< [ [ \G(x,y,t,s)\\h(y,s)\dyds. (13) 

do Jq 

Also note that G satisfies the inequality (see, e.g., [77] 1 

\G(x,y,t,s)\ < x,y £ Q, 0 < s < t, (14) 

where C\ > 0 does not depend on x, y, t or s. 

Similarly to the proof of Theorem 12.41 for every s < t the integral of \h(y, s)| 
over Q is bounded by ||&i — &2||c[o,t] an< i hence by ||ai — O2||c[o,t] an d hence by 
||ui — ^ 2 IIc(q7) • Combining this with (fT51) and (THl) . and taking the supremum 
over (, x , t) £ Qt we get 


IMIc«t) — C>VTMc iQT) , 

where C 2 > 0 does not depend on T. Thus w = 0 for T small enough. A passage 
to arbitrary T is standard. 
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Theorem 2.6 (See |11L Theorem 2.9]) Let u £ be a transverse 

solution of problem © ©• If \[P — p n \\ w a-2/«, — )• 0 and \b n — b\ — > 0 as n —> 

W q,N WJ 

oo, then for sufficiently large n, problem ©~([8]l has a solution u n £ Wg' 1 (QT) 
with initial data p n and initial configuration fon defined via b n . Furthermore, 
\\u n (Qt) ^ ^ as n t oo. 

The crux of the proof is showing that for sufficiently large n, all the solutions 
exist on the same time interval [0,T]. To this end we note that we have in fact 
given an explicit construction of T, and that this T depends on b , /i^gy 

and if p[b) = /3, also on p'(b). Hence for p n and b n close enough to p and b in 
their respective norms, the same T can be used. 


2.3 Case n > 2 

For the case n > 2 a notion of transversality has been studied in a model prob¬ 
lem. For clarity we will define transversality for the case where the threshold j3 
is adjoined to the free boundary between Q i and Q-i, and a is not. The gen¬ 
eral case is treated similarly. In what follows, let int(H) denote the topological 
interior of a subset A C Q, and let {p = a} be defined similarly to {u = a} 
but taking x £ Q instead of (x,t) £ Qt- In [ 22 ] the existence and uniqueness 
of solutions were studied for initial data transverse in the following sense (see 
Figure HJi and0]D, and recall that Qj is given by (©). 

Definition 2.7 We say the function p is transverse with inspect to £o if the 
following hold: 

(1) Q i and Q-i are measurable, dQ -1 C Q, dQ\ = <9Q-i U dQ, and dQ-i 
has zero Lebesgue measure. 

(2) p(x) < ft for x £ int(Qi) U dQ. 

(3) p{x) > a for x £ Q-\. 

(4) If x £ (u = /3} fl dQ-i, then there is a neighbourhood A of x, a set 
A! C K" -1 , a k > 0, and a map if such that 

(4- a) if is a composition of a translation and a rotation and 

if (A) = A! x [— k, «], if(x) = (0,0). 


(4-b) There is a continuous function b : A! —> [— k, «] such that the con¬ 
figuration function fo o ifin if (A) (which we denote by fo{y' ,Vn), 
y' £ A') is given by 


fo{y',y n ) 


-1 if - k < y n < b{y’), 
1 */ b(y') < y n < k. 


(4-c) p o if which we write as p(y',y n ), satisfies p y „(0,0) < 0. 


10 




Figure 4: An example of a) the sets Q±i, b) transverse data, and c) non- 
transverse data. 


We observe that in Subsection 12.21 the boundary between Q\ and Q-\ was a 
single point b. But when n > 2, this boundary is assumed to have the structure 
of a continuous codimension 1 submanifold in a neighborhood of a point on the 
free boundary where <p takes a threshold value. Also note that for n > 2 non- 
transversality can be caused by the geometry of <9Q-i in addition to the possible 
degeneracy of SAp (see Figureand Subsection 12.41 for further discussion). 

Theorem 2.8 (see (221 . Theorems 3.18 and 3.19]) Assume that n > 2 and 
T 6 < N^ q (Q) is transverse with respect to £o- Then there is a T > 0 such 
that any solution u £ W^ ,:l (Qt) to problem (|6|)- (j8jl is transverse and there is 
at least one such solution. Moreover, if for some T' > 0, ui and U 2 are two 
transverse solutions to problem © © on Qt> , then u\ = U 2 - 

The main ideas of the proof are similar to those for the case n = 1 . Since 
•), £o(y', •)) is transverse in the Id sense for every y' £ A', one can prove 
continuity of a map 1Z that now maps functions u o £ C x (Qt ) (A < 7) to 
solutions TZ(uo) := u of problem © © with the right hand side %(£o,«o)- 
Estimate m implies that u £ C^^Qt), and the compactness of the embedding 
C 7 (Qt) C C x (Qt) and the Schauder fixed point theorem together imply that 
1Z has a fixed point in C 7 (Qt)- 

2.4 Generalizations and Open Problems 

Let us list some generalizations for the case n = 1. 

Change of topology. Suppose u(x, t) becomes non-transverse at some 
time T in the sense of Definition 12.21 Then one of two possibilities arise. Either 
u(x, T) has touched a threshold with zero spatial derivative at some point in 
(0,1), or this is not the case but lim t ^ T b{t) = 1. In the latter case, one 
can continue the solution, and it remains unique, by redefining the problem 
effectively without hysteresis ED Theorem 2.18]. We say that the topology 
of the hysteresis has changed at time T, in the sense that £ transitions from 
piecewise constant to uniformly constant. 
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Continuous dependence on initial data. If u is a solution such that the 
topology has changed for some t\ < T, then u need not continuously depend 
on the initial data since a sequence of approximating solutions u n may become 
non-transverse at moments r n with r n < t\ and limn-^ r„ = t\ (the dashed 
line in Figure [5]). But if we also assume that each u n is a transverse solution, 
then solutions do depend continuously on their initial data. 


£(avr 0 ) £(x ,r n ) £0Mi) 


1 

r- 1 

1 6(to) ] 

'b{T n ) 1 

1 

T 

* 

t 1 * 

1 

1 

-1 





Figure 5: A solution u (drawn as solid lines in the lower picture) and its con¬ 
figuration £ (the upper picture) that remain transverse as a discontinuity of £ 
disappears at time t\. The dashed line in the lower picture is a series of non- 
transverse approximations u n that become non-transverse at moments r„ with 
r n < ti and lirn^oo r n = t\. 


Finite number of discontinuities. The results in Subsection 12.21 remain 
valid if the hysteresis topology is defined by finitely many discontinuity points. 
The hysteresis changing topology in the sense we described for one point of 
discontinuity corresponds to these points merging together in the general case 
(see Figure [6]). 

General nonlinearity. The results in this section also hold for the more 
general problem 0 ©■ First one must assume that / is locally Lipschitz and 
dissipative (see m Condition 2.11]). With such an /, and if Hi and if_i are 
locally Holder continuous, then transverse solutions exist and can be continuted 
up to a maximal interval of transverse existence. If one additionally assumes 
that transverse solutions are unique, they can be shown to continuously depend 
on their initial data. To prove the uniqueness of solutions the authors of [201122] 
make the stronger assumption on Hi and H_ i, namely that 


\Hi(ui) 


Hi(u 2 )\ < 


M 

(/? - my + (/3 - u 2 y 


\ui 


U 2 \, 


for ui,u 2 in a left neighborhood of [3, with M > 0 and a £ (0,1), plus an 
analogous inequality for _H_i and a right neighborhood of a. This condition 
covers the case where Hi and ff_i are the stable branches in the slow-fast 
approximation as in Figure [2] (see the appendix of 1201 for further discussion). 
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£(Mo) £(M») f(Mi) 



Figure 6: Discontinuities merging as t —> t\. 


Systems of equations. In [51:, Theorem 2.1], the results of Subsection l2.2l 
were generalized to systems of equations of the type in problem © ©• It was 
also shown therein that problem QJ ([3]) can be coupled to ordinary differential 
equations to cover the Hoppensteadt-Jager model from mm- 

Let us conclude this subsection by discussing an open problem. 

Open problem. In Figure [T]\ one can see that for every y' ^ 0, 
•), £o(y\ •)) ^ transverse in the Id sense (with two discontinuties), but 
since the free boundary cannot be represented as a graph with codomain y n at 
the point y' = 0, this initial data is not transverse. Whether Definition 12.71 can 
be generalized to include such cases is the subject of future work, and at this 
stage the authors strongly suspect that item [4] of Definition 12.71 can be replaced 
by the following statement: if x £ {u = /3} fl dQ-i, then S7ip{x) / 0. In other 
words, the assumption that the free boundary is a graph is not necessary, and 
hence Figure would also be transverse. This question is intimately linked to 
the topology of the free boundary. Whether solutions can be continued to a 
maximal interval of existence and how to pose continuous dependence of initial 
data is unclear for the quite general conditions on Q-\ and Q\ in Definition 12.71 
These questions also apply to the case where n = 1 and has infinitely many 
discontinuities. 


3 Regularity of Strong Solutions 

To begin with let us discuss what we mean by regularity of solutions in this 
context. First observe that we cannot expect a classical solution since TL has a 
jump discontinuity. Therefore the “optimal” regularity we expect is W^ 1 . In 
this section we obtain W^ 1 “locally”, for points ( x , t) € Qt outside of the static 
part of the free boundary. We will also assume the following condition: 

Condition 3.1 Hi(u) = 1 and H_i(u) = — 1. 
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Let us introduce the notation Q^ 1 := {(a;, t) | £{x,t) = ±1} and observe that u 
is smooth on the interior of QIf 1 . 

The free boundary is defined as the set V := dQlp (~l OQf . Moreover, we 
define T a := {u = a} D T and := {u = (3} fl T. Note that both T Q and 
have zero Lebesgue measure whenever u is a solution of problem (0) (151) . This 
follows from the fact that ut — A u = 0 a.e. on T Q U T^ and Condition 13.11 (see 
Alt’s argument in the introduction and mr 

The estimates we obtain will depend critically on the static part of the free 
boundary T„ := r\(T Q ur^). If ( x, t) € F„, then u(x, t) ^ a,j3 and by continuity 
of u, u(x, t ± r) ^ a, p for r sufficiently small. This means £(ir, t ± r) = £(x, t) 
and so if we draw the Taxis vertically as in Figure 0 T„ looks like a vertical 
strip. 



Figure 7: A possible scenario where T / F a U and r„ appears. White and 
grey indicate the regions Q? and Qf 1 respectively. 

Next we recall the definition of a parabolic cylinder 

P r (x°,t°) := {ieK"j ||x° - a:|| R » < r} x (t° - r 2 , t° + r 2 ), r> 0. 

We define the parabolic distance between (x°,t°) and a set A C Qt as 

dist p ((x 0 , t°), A) := sup{r > 0 | P r (x°,t°) n {t < £ 0 } fl A = 0}. 

This is all the notation we need to state the main result of [23]. 

Theorem 3.2 (see |23l Theorem 2.3]) We assume that n > 1 and u is a 
solution of problem (0) (|8|). Then 

n 

\ut(x,t)\ + ^2 \u XiXj {x,t)\ < C(p v ,pb,M), a.e. (x,t) € 

*,j= 1 
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where C depends on p v := dist p ((x,t),T v ), p b := dist p ({x, t), d'Qr U (Q x {0})), 
and M := sup^ t)eQT \u(x, t)\. 

To explain the main ideas in the proof we define some further notation. Let 
r° = r Q fl {Vu = 0} and T* = r a \r°, with rjg and defined similarly. 
Furthermore, define r° = U and T* = F* U 

The crucial point in the proof is the quadratic growth estimate 

sup \u - /3| < Ci (p v , P b,M)r 2 for r < minjp,,^;,} , (15) 

P r (x,t) 

and (a;,f) £ T° (the estimate on T° is similar). The main tool for showing the 
quadratic bound (1151) is the local rescaled version of the Caffarelli monotonicity 
formula, see [23U281f29] . 

Furthermore, the quadratic growth estimate (fT51) implies the corresponding 
linear bound for |Vzt| 

sup |Vu| < C- 2 {p Vl Pb, M)r for all r < min{jO„,p b } , (16) 

P r (x,t ) 

with (x,t) € r°. The dependence of C\ and C 2 on the distance p v in CSD and 
m arises due to the monotonicity formula. Near neither the local rescaled 
version of Caffarelli’s monotonicity formula nor its generalizations (such as the 
almost monotonicity formula) are applicable to the positive and negative parts 
of the spatial directional derivatives D e u , with e £ R". 

Besides estimates m and m, one also needs information about the be¬ 
haviour of lit near T*. Although ut may have jumps across the free boundary, 
one can show that ut is a continuous function in a neighborhood of (a;, t) £ 
T* \ r„. In addition, the monotonicity of the jumps of 'H(^opu) in the t- 
direction provides one-sided estimates of ut near r a and Combining these 
results with the observation that u t < 0 on T* \ and u t > 0 on F)j \ gives 

sup \u t \ < C 3 (p bl M). (17) 

r*\r„ 

Inequalities m o allow one to apply methods from the theory of free bound¬ 
ary problems (see, e.g., [11 s TQ3j and estimate \u t (x,t)\ and \u XiXj (x,t)\ for a.e. 

(x,t) £ Qt\ r„. 


4 Non-Transverse Initial Data 

4.1 Setting of a Problem 

In this section we summarize the recent work [24] , where the nontransverse case 
is analyzed for x £ R, and indicate directions for further research. We will 
be interested in the behavior of solutions near one of the thresholds, say /3. 
Therefore, we set a = —00 and /3 = 0 (see FigureH]) and assume that the initial 
data satisfy ip(x) = — cx 2 + o(x 2 ) in a small neighborhood of the origin, ip(x) < 0 
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H(to,u) 


A 

hi 


P = 0 


u 

► 


Figure 8: Hysteresis with thresholds a = — oo and ft = 0. 


everywhere outside of the origin, £o(x) = —1 for x = 0, and £o{x) = 1 for x yf 0. 
In particular, we assume c > 0. In this situation, the theorems in Section 12.21 
are not applicable. Hence, to understand the dynamics of the solution near the 
origin, we approximate the continuous equation © by its spatial discretization 
and the initial data by the discrete quadratic function. Namely, we choose a grid 
step e > 0, set u e n (t) := u(sn,t), n G Z, and consider the system of infinitely 
many ordinary differential equations with hysteresis 


du £ n 

dt 


X 'n+1 


- 2< + ul 


+ hi{u c n ), 


t > 0, n G Z, 


(18) 


supplemented by the nontransverse (quadratic) initial data 


u e n (0) = — c(en) 2 , n G Z. 


(19) 


Here we do not explicitly indicate the dependence of W on (o, assuming that 
TL(iin)(t) = hi if Un(s) < 0 for all s G [0,i] and H(u^){t) = h-i otherwise. As 
before, we assume that h_i < 0 < hi. 

Due to [24], Theorem 2.5], problem (fT8l) . (fl9l) admits a unique solution in the 
class of functions satisfying 

sup K(s)| < Ae B \ n \, n € Z, t> 0, 
se[o,t] 


with some A = A(t, e) > 0 and B = B(t, e) G R. Thus, we are now in a position 
to discuss the dynamics of solutions for each fixed grid step e and analyze the 
limit £ —^ 0. 

First, we observe that e in (fT8]l . (fT9l) can be scaled out. Indeed, setting 

u n (t) := £~ 2 u e n {E 2 t) (20) 


reduces problem (TT51) . (fTHl) to the equivalent one 


du n 

dt 


— Un+l + U n — i 


B-iUn), 


t > 0, n G Z, 


( 21 ) 


,(0) = -c 


n G Z. 
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Using the comparison principle, it is easy to see that if hi < 2c, then u n (t) < 
0 for all n £ Z and t > 0 and, therefore, no switchings happen for t > 0. Let us 
assume that 

h -1 < 0 < 2c < hi. (22) 

It is easy to show that u n (t) < 0 for all n £ Z and t > 0. However, some nodes 
can now reach the threshold 0 = 0 and switch the hysteresis. The main question 
is which nodes do this and according to which law. 


4.2 Numerical Observations 

The following pattern formation behavior is indicated by numerics (see Fig¬ 
ure 0). As time goes on, the spatial profile of u n [t) forms two symmetric hills 



Figure 9: Upper graphs represent spatial profiles of the hysteresis H(m„) and 
lower graphs the spatial profiles of the solution u n . a) Nontransverse initial 
data, b) Spatial profiles at a moment t > 0 for h—i = 0. c) Spatial profiles at a 
moment t > 0 for h-i = —hi < 0. 


propagating away from the origin. At the same time, the whole spatial profile 
oscillates up and down (never exceeding the threshold 0 = 0) and touches the 
threshold 0 = 0 in such a way that 


lim 

j-* °o 


NnsU) 

N s (j) 


\h-i\ 

hi 


(23) 


where N s (j) and N ns (j) are integers denoting the number of nodes in the set 
{mo>w±i, ... ,u±j} that switch and do not switch, respectively, on the time in¬ 
terval [0,oo). In |24] . such a spatio-temporal pattern was called rattling. 

A more specific pattern occurs if \h-i\/hi = p n s/ps, where p s and p ns are co¬ 
prime integers. In this case, for any j large enough, the set {wj+i, • • •, Uj+p s +p ns } 
contains exactly p s nodes that switch and p ns nodes that do not switch on the 
time interval [0,oo). 

If a node u n switches on the time interval [0, oo), then we denote its switching 
moment by t n ; otherwise, set t n := oo. In particular, finite values of t n char¬ 
acterize the propagation velocity of the two hills mentioned above. Numerics 
indicates that, for the nodes where t n is finite, we have 


t n = an 


2 



if i = 0, 

if i < 0, 


as n —> oo, 


(24) 
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and 

\uk+i(t) - Uk(t)\ < b, \k\ < n, t > t n , n = 0,1,2,..., (25) 

where a, b > 0 do not depend on k and n. In particular, (El) and (El) mean that 
the hills propagate with velocity of order f -1 / 2 , while the cavity between the 
hills has a bounded steepness, which distinguishes the observed phenomenon 
from the “classical” traveling wave situation. 

4.3 Rigorous Result 

The recent work [241 provides a rigorous analysis of the rattling in the case 
h _i = 0, where, according to (El) , all the nodes are supposed to switch at time 
moments satisfying 

t n = an 2 + q n , \q n \ < Eyfn , (26) 

where E > 0 does not depend on n £ Z. In |24j . the authors found the coefficient 
a and proved that if finitely many nodes u n (t ), n = 0, ±1, • ■ ■ ± no, switch at 
time moments t n satisfying (1261) . then all the nodes u n (t), n £ Z, switch at 
time moments t n satisfying (Ell (see the rigorous statement below). One of the 
main tools in the analysis is the so-called discrete Green function y n (t) that is 
a solution of the problem 


( Vo = &y 0 + 1, t > 0, 

< y n = A y n , t > 0, n/0, (27) 

[yn(0) = 0, n G Z. 

The important property of the discrete Green function is the following asymp¬ 
totics proved in [30] : 

y n (t) = + O ^—0 as t —>• oo, (28) 

where 

OO 

f{x) := 2x J y~ 2 h(y) dy, h(x) := — e~^, (29) 

X 

and O(-) does not depend on n € Z. 

Now if we (inductively) assume that the nodes u o, u± i,... u±( n _ i) switched 
at the moments satisfying while no other nodes switched on the time 

interval [0,£ n _i], then the dynamics of the node u n (t) for t > t n _\ (and until 
the next switching in the system occurs) is given by 

n— 1 

U n (t) = -cn 2 + (hi - 2c)t - hi ^ y n - k {t-t k ). (30) 

k=—(n— 1) 
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At the (potential) switching moment t n = an 2 + q n , the relations tk = ak 2 + q & 
(jfcj < n— 1), equality (l30l) . the Taylor formula, and asymptotics (l28l) yield 

n— 1 

0 = —cn 2 + (hi — 2 c)an 2 — hi ^ y n -k (a(n 2 — k 2 )) + l.o.t. 

k=—(n— 1 ) 


n— 1 


= —cn 2 + (hi — 2 c)an 2 — hi ^ \Ja(n 2 — k 2 )f j 

k=—(n— 1 ) 


n — k 


y/a(n 2 — k 2 ) 


(31) 


+ l.o.t. 


= (—c + (hi — 2 c)a — hiR n (a )) n 2 + l.o.t., 


where 


n ~ 1 i ,_ / 

Rn(a):= ~\/a(l- (fc/n) 2 )/ 

Z ' n \ 

k=-(n- 1 ) \ 


1 — A:/n 


\/a(l - ( k/n ) 2 ) 


and “l.o.t.” stands for lower order terms that we do not explicitly specify here. 
Note that R n (a ) is the Riemann sum for the integral 


!,(“) := / ^( 1 -» 2 )/( y „ 1 (1 ( C ia) ) < 32 ) 

Therefore, equality m can be rewritten as 

0 = (—c + (hi — 2 c)a — hilf(a )) n 2 + l.o.t. (33) 

It is proved in |24] that there exists a unique a > 0 for which the coefficient at 
n 2 in (1331) vanishes. The most difficult part is to analyze the lower order terms 
in (l33l) that involve: 

(1) the remainders qo, q±i ,..., q n from (l26l) arising from (l30l) via the applica¬ 
tion of the Taylor formula, 

(2) the remainder in the asymptotic (l28l) for the discrete Green function y n (t ), 

(3) the remainders arising from approximating the integral If (a) by the Rie¬ 
mann sum R n (a). 

In particular, one has to prove that if \qj\ < Ey/\j\ for j = 0, ±1,...,±(n — 1), 
then the lower order terms vanish for a specified above and \q n \ < Eyf\n\. This 
allows one to continue the inductive scheme and (after an appropriate analysis 
of the nodes u±( n+ i)(t),a±( n+2 )(t), ... for t £ [t n -i,t n ]) complete the proof. 
The rigourous formulation of the main result in [2!] is as follows. 

Theorem 4.1 (see j24|, Theorem 3.2]) Assume that (l22l) holds and that 
h_i = 0. Let a = a(hi/c ) > 0 be a (unique) root of the equation 


— c + (hi — 2 c)a — hilf(a) = 0 


(34) 
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with If (a) given by (1351) . Then there exists a constant Eo = Eo(hi,c,a ) > 0 
and a function no = no(E) = no(E,hi,c,a) (both explicitly constructed ) with 
the following property. If 

finitely many nodes uo(t),ui(t ),..., u no (t) switch at moments t n 

(35) 

satisfying (1261) with the above a and some E > Eq, 


then each node u n (t), n € Z, switches; moreover, the switching occurs at a time 
moment t n satisfying m with a and E as in (ESD- 

We note that the explicit formula (l30l) for the solution u n (t) allows one to 
verify the fulfillment of finitely many assumptions (1351) numerically with an ar¬ 
bitrary accuracy for any given values of hi and c. The graphs in Figure [TU] taken 
from j53] represent the values of a, E, and no(E) that fulfill assumption (1551) 
for c = 1/2 and hi = 1.1,1.2,1.3,..., 2.5. 




Figure 10: Dependence on hi of the values of a, E, and no(E) that fulfill assump¬ 
tions (l35l) for c = 1/2. a) The values of a are found as roots of (l34l) . b), c) The 
values of E and no{E) are calculated for discrete values hi = 1.1,1.2,1.3,..., 2.5. 


4.4 Open Problems 

To conclude this section, we indicate several directions of further research in the 
nontransverse case. 

Case h- 1 < 0. In this case, one has to additionally prove a specific switching 
pattern (l23l) . We expect that the tools developed in [24] will work for rational 
hi/h-i. The irrational case appears to be a much more difficult problem. 

Multi-dimensional case. Numerics indicates that the behavior analogous 
to (1231) occurs in higher spatial dimensions for different kinds of approximating 
grids. Figure |TT] illustrates the switching pattern for a two-dimensional analog 
of problem (121[) , where the Laplacian is discretized on the square and triangular 
lattices, respectively. 

Limit e —> 0. We introduce the function 

u e (x, t) := u e n (t), x G [en — e/2, en + e/2), n € Z, 

(which is piecewise constant in x for every fixed t ). Making the transformation 
inverse to m and assuming (HU) and (l24l) . we can deduce that, as e —> 0, 
the function u E {x,t) approximates a smooth function u(x,t ), which satisfies 
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Figure 11: A snapshot for a time moment t > 0 of a two-dimensional spatial 
profile of hysteresis taking values hi > 2c > 0 and h_i = —hi < 0. The 
nontrasverse initial data is given by ip(x) = — c{x\ + x 2 )- Grey (black) squares 
or hexagons correspond to the nodes that have (not) switched on the time 
interval [0, £]. a) Discretization on the square lattice, b) Discretization on the 
triangular lattice. 


u{x,t) = 0 for x £ (—y/t/a, y/t/a). In other words, u[x,t) sticks to the thresh¬ 
old line /? = 0 on the expanding interval x £ (~y/ t/a , y/t/a). 

Similarly to u e (x,t ), we consider the function 

H £ (x,t) :=H(u%)(t), x £ [en — e/2, en + e/2), n £ Z, 

which is supposed to approximate the hysteresis H(u)(x,t ) in We see that 
the spatial profile of H £ (x,t) for x £ {—y/t/a , y/t/a) is a step-like function 
taking values hi and h_i on alternating intervals of length of order e. Hence, it 
has no pointwise limit as e —> 0, but converges in a weak sense to the function 
H(x,t) given by H(x,t) = 0 for x £ (— i/t/a, y/t/a) and H(x,t) = hi for x ^ 
{—yjt/a, y/t/a). We emphasize that H(x,t ) does not depend on h- 1 (because 
a does not). On the other hand, if h- 1 < 0, the hysteresis operator 'H(u){x,t) 
in m cannot take value 0 by definition, which clarifies the essential difficulty 
with the well-posedness of the original problem (jGj) in the nontransverse case. 
To overcome the non-wellposedness, one need to allow the intermediate value 0 
for the hysteresis operator, cf. the discussion of modified hysteresis operators 
due to Visintin and Alt in the introduction. A rigorous analysis of the limit 
6 —^ 0 is an open problem, which may lead to a unique “physical” choice of 
an appropriate element in the multi-valued Visintin’s hysteresis Hvis(£o> M ) i n 
Definition 11.21 

Rattling in slow-fast systems. One may think that the rattling occurs 
exclusively due to the discontinuous nature or hysteresis. This is not quite the 
case. Consider an equation of type (0 with the hysteresis %(£o,m) replaced by 
the solution v of a bistable ordinary differential equation of type (J4]) , e.g., 

ut = u X x + v, Sv t = g(u,v). (36) 

Numerical solution of system (13611 with a nontransverse initial data u(x, 0) = 
—cx 2 + o(x 2 ) and v(x, 0 ) = Hi(/3) near the origin reveals a behavior analogous 
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to that for a spatially discrete system (see Figure H2l). As the spatial profile of 
u(x,t) touches the threshold /? at some point xo, the spatial profile of v{x,t) 
forms a peak-like transition layer around Xq that rapidly converges to a plateau. 
Thus, as time goes on, the spatial profile of v(x,t) converges to a step-like 
function taking values Hffj3) and FZ_i(/3) on alternating intervals, whose length 
tends to zero as 8 —> 0. A rigorous analysis of the limit S —> 0 is an open 
problem. 

tends to 0 

as 6 —>• 0 



Figure 12: Lower and upper graphs are spatial profiles of the solution u{x,t) 
and v(x,t), respectively, for problem (l36l) with initial data m|*=o = — ca,’ 2 + o(ar), 

dt=o = 
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